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Abstract. We present the results of a new set of hydro- 
dynamical simulations of hot gas flows in model ellipti- 
cal galaxies with the following characteristics: the spatial 
luminosity distribution approaches a power law form at 
small radii, in accordance with the results of recent ground 
based observations, and with the Hubble Space Telescope; 
the dark matter has a peaked profile too, as indicated 
by high resolution numerical simulations of dissipation- 
less collapse; the dark to luminous mass ratio spans a 
large range of values, including low values found by op- 
tical studies confined to within two effective radii; finally, 
the type la supernova rate is that given by the latest esti- 
mates of optical surveys, or zero, as suggested by the iron 
abundances recently measured in the hot gas. 

We find that the resulting gas flows are strongly decou- 
pled: an inflow develops in the central region of the galax- 
ies, while the external parts are still degassing, i.e., the 
flows are mostly partial winds. This behavior can be ex- 
plained in terms of the local energy balance of the hot gas. 
A large spread in the X-ray luminosity Lx at fixed optical 
luminosity Lb can be produced as in previous simulations 
that used King models plus massive quasi-isothermal dark 
halos, and higher supernova rates; the key factor causing 
large Lx variations is now the size of the central inflow 
region. The highest Lx observed correspond to global in- 
flows. Finally, non negligible amounts of cold gas can be 
produced by the partial winds; this could be an explana- 
tion for the possible discovery of cold matter at the center 
of elliptical galaxies, an alternative to the presence of a 
steady state cooling flow. 

Key words: Galaxies: elliptical and lenticular, cD - 
Galaxies: ISM - Galaxies: kinematics and dynamics - X- 
rays: galaxies 
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1. Introduction 

X-ray observations, beginning with the Einstein Observa- 
tory, have demonstrated that normal early-type galaxies 
are X-ray emitters, with 0.2-4 keV luminosities ranging 
from ~ 10 40 to - 10 43 erg s" 1 ( |Fabbiano 1989] Fab- 
biano, Kim, & Trinchieri 1992). The X-ray luminosity 
Lx is found to correlate with the blue luminosity Lb 
(L x oc L B 0±0 ' 2 ), but there is a large scatter of roughly two 
orders of magnitude in Lx at any fixed Lb > 3 x 1O 1O L0 
(Fig. 1). The observed X-ray spectra of galaxies with high 
Lx/Lb ratios are consistent with thermal emission from 
hot, optically thin gas, while those of low Lx/Lb ob- 
jects can be mostly accounted for by emission from stellar 
sources (Kim, Fabbiano, & Trinchieri 1992). 

The scatter in the Lx — Lb diagram has been rec- 
ognized as the most striking feature of the X-ray proper- 
ties of early- type galaxies. Using new apparent magnitudes 
and fundamental plane distances, it was shown that this 
scatter is reduced by 20%, but not eliminated (Donnelly, 
Faber & O'Connell 1990; see Fig. lb). This result is based 
on the old estimate of the X-ray fluxes (that of Canizares, 
Fabbiano, & Trinchieri 1987), and on just half of the final 
sample of X-ray galaxies produced by Fabbiano et al. 1992 
(see Eskridge, Fabbiano & Kim 1995 for a detailed com- 
parison of the statistical results by Donnelly et al. with 
those obtained using the whole sample). One can argue 
that the large dispersion in Lx is definitively not the result 
of distance errors on the basis of the fact that a scatter of 
the same size as in Fig. la is present even in the distance- 
independent diagram of Lx/Lb versus the central stellar 
velocity dispersion (e.g., Eskridge et al. 1995). 

Many theoretical models were developed to explain 
the findings above, including numerical simulations of 
the behavior of gas flows fed by stellar mass loss and 
heated by type la supernovae (SNIa). Steady state cool- 
ing flow models were in vestigated first (Nulsen et al. 1984, 
Barazin fc White 1987 ^988), and it was found that these 
can only reproduce X-ray bright galaxies. Evolutionary 
models with a SNIa rate approximately constant with time 
have been carried out by Mathews & Loewenstein (1986), 
Loewenstein & Mathews (1987), and David, Forman, & 
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Jones (1990, 1991). After a very brief initial wind phase, 
driven by the explosion of type II supernovae, the result- 
ing flow evolution goes from a global inflow to a wind, 
which is experienced only by the smallest galaxies by the 
present time. David et al. (1991) conclude that all galax- 
ies above Lb — 3 x 10 10 L© host a cooling flow. So, as 
in the steady state cooling flow scenario, the scatter in 
Lx at fixed Lb has to be explained by a combination of 
environmental differences ( White fc Sarazin 1991 ) and by 
large variations from galaxy to galaxy in the stellar mass 
loss rate per unit Lb, the efficiency of thermalization of 
the stellar mass loss, and that of thermal instabilities in 
the hot gas. An alternative way of explaining the scatter in 
the Lx — Lb diagram is given by the evolutionary scenario 
of D'Ercole et al. (1989), and Ciotti et al. (1991, CDPR). 
Assuming that the SNIa explosion rate is declining with 
time slightly faster than the rate with which mass is lost by 
the stars, in the beginning the energy released by SNIa's 
can drive the gas out of the galaxies through a supersonic 
wind. As the collective SNIa energy input decreases, a 
subsonic outflow takes place, which gradually slows until 
a central cooling catastrophe leads to the onset of an in- 
flow. At fixed Lb, any of the three phases wind, outflow or 
inflow can be found at the present epoch, depending only 
on the various depths and shapes of the potential well of 
the galaxies. In this way both the large scatter in Lx and 
the trend in the spectral properties are accounted for at 
the same time: in the X-ray bright galaxies the soft X-ray 
emitting gas dominates the emission, being in the inflow 
phase, that resembles a cooling flow; in the X-ray faint 
galaxies the hard stellar emission dominates, these being 
in the wind phase. 

Recent ASCA observations seems to indicate that the 
SNIa's activity is suppressed in early-type galaxies, which 
implies that the CDPR scenario is essentially ruled out. 
Since this issue is far from closed (§1.3), it is still worth- 
wile to explore the effects of SNIa's on the flows, using 
the updated rate given by recent optical surveys; this is 
lower than that adopted by CDPR, which was 0.67 the 
standard one estimated by Tammann (1982). This aspect, 
together with the need for changes in other galaxy proper- 
ties crucial for the evolution of hot gas flows, are discussed 
below. 

1.1. New stellar density profiles 

Previous numerical simulations of hot gas flows used the 
King (1972) stellar density distribution mainly for com- 
putational ease. This distribution has an inner region of 
constant density (the so called "core" ) of the order of a 
few hundred parsecs, which keeps the time step of the 
numerical simulation reasonably small (a stellar density 
increase implies a reduction of the characteristic hydrody- 
namical time step) . Another advantage of density profiles 
with a constant density core is that the central regions 
can be resolved using a larger grid size; this again allows 



Fig. 1. The Lx — Lb diagram of early- type galaxies with X-ray 
emission detected by the Einstein satellite; X-ray fluxes are 
from the final catalog of Fabbiano et al. (1992). Apparent B 
magnitudes and distances are from Fabbiano et al. (1992) in 
Fig. la, and from Donnelly et al. (1990) in Fig. lb (see Section 
1). The dashed line gives an estimate for the stellar source 
contribution to Lx (from Kim et al. 1992). Also shown with 
triangles are the positions of the models calculated here (see 
Section 3.5), to which the stellar source contribution has been 
added. 



a larger time step (see CDPR for a more quantitative dis- 
cussion). It is well known though that a much better de- 
scription of the surface brightness profiles of ellipticals is 
given by the de Vaucouleurs (1948) law. A very good fit 
to this law is given by the Jaffe (1983) and by the Hern- 
quist (1990) distributions, that have the advantage that 
all their dynamical properties can be expressed analyti- 
cally. These distributions belong to the family of the so 
called 7-models, that has been widely explored recently 
( pchncn 199^ ; [Tremaine et al. 1994[ ): 



p{r) 



M— 



7) 



47r ri(r c +r) 4 -i' 



(1) 



The density profile of the Hernquist law has 7 = 1, 
while that of the Jaffe law has 7 = 2. Inside r c the density 
of 7-models increases as r -7 , a significant difference with 
respect to King models: over a few hundreds of parsecs at 
the center the 7-models are power laws. The existence of 
cores of constant surface brightness has been definitively 
ruled out by ground based observations (Lauer 1985, Ko- 
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rmendy 1985), and recently by the Hubble Space Tele- 
scope (Jaffe et al. 1994; Lauer et al. 1995), that has shown 
how the central surface brightness profile is described by 
a power law as far in as can be observed, i.e., ~ 10 pc in 
Virgo. From the point of view of a very accurate photom- 
etry, the surface brightness law implied by eq. (1) cannot 
reproduce well both the envelope and the very center of 
all the ellipticals observed by HST. The Jaffe law, though, 
gives a general description of ellipticals accurate enough 
for the treatment of hot gas flows, on the scales that are 
relevant for the problem (from a few tens to several thou- 
sands of parsecs). 



retical arguments favor a peaked profile (Ciotti & Pel- 
legrini 1992; Evans & Collett 1997), and high resolu- 
tion numerical simulations of dissipationlcss collapse pro- 
duce a density distribution with 7 ~ 1 near the cen- 



ter (Dubinski & Carlberg 1991; Navarro, Frenk, & White 



1996; White 1996; Fukushige & Makino 1997, and refer- 



ences therein). Previous works studying the hot gas flow 
evolution always used quasi isothermal haloes at least nine 
times more massive than the luminous component. We are 
motivated now to explore even the effects produced by 
dark haloes not as massive as supposed before, and not 
quasi-isothermal . 



1.2. New dark matter estimates and distributions 

Attempts to estimate the amount of nonluminous mass in 
elliptical galaxies have been made recently through exten- 
sive searches for dynamical evidence of dark matter, either 
with systematic ob servations of ionize d gas probing the 
gravitational field ( Pizzclla et al. 1997 ) or with measure- 
ments of s tellar velocity disp ersions profiles out to large 
radii (e.g., [Bertin et al. 1994| Carollo et al. 1995). These 



optical studies are confined to within one or two effective 
radii i? c , and typically find that luminous matter domi- 
nates the mass distribution inside ~ R c , while dark matter 
begins to be dynamically important at 2-3 R c . In partic- 
ular, for a sample of X-ray emitting galaxies, it has been 
found that dark matter halos are not much more massive 
than the luminous component, with the dark-to-luminous 
mass ratio Mh/M* rs 1 — 6; the value of Mh/Af* ~ 2 



is most common (Saglia et al. 1993). X-ray emission from 



hot gas provides a great potential for mapping the mass 
of ellipticals to larger distances (e.g., Fabian et al. 1986). 
The standard method employed derives from the equa- 
tion of hydrostatic equilibrium, and requires the knowl- 
edge of the gas temperature profile. Attempts to apply 
this technique to the Einstein data yielded a much larger 
component of dark matter than found from optical data, 
but these results are very uncertain because temperature 
profil es are poorly de termined (Forman, Jones, & Tucker 
1985; Fabbiano 1989 ). Using superior X-ray data provided 
by ROSAT and ASCA, Buote & Canizares (1997) proved 
that mass does not follow the optical light, out to many 
i? c , in NGC720 and NGC1332. Adopting plausible gas and 
mass models, they find M h /M* > 3 for NGC1332, and 
M h /M* > 7 in NGC720 at 90% confidence; M h prevails 
exterior to R c . Similarly Mushotzky et al. (1994) derive 
M h /M* ~ 8 within 8_R G for NGC4636, and analogous re- 
sults have been obtained from ROSAT data of NGC507 
and NGC499 by Kim & Fabbiano (1995). AXAFvnR have 
the combined spatial and spectral resolution to measure 
accurately the presence of different spectral components, 
their relative flux and their spatial distribution, and to 
produce more accurate mass distribution from X-rays. 

The radial density distribution of the dark haloes of 
ellipticals is not well constrained by observations; theo- 



1.3. New Type la supernova rates 

Nearby SNIa rates in early-type galaxies have been 
carefully reanalyzed recently, and this important in- 
gredient of the simulations of hot gas flows has 
been revised. From optical surveys it was esti- 
mated to be 0.88 h 2 SNu (fTammann 1982|), an d then 
0.98 h 2 SNu ( |van den Bergh fc Tammann 199l[ ). Most 
recent estimates agree on lower values: 0.25 h 2 SNu 



( |van den Bergh fc McClure 1994j ) and 0.24 h 2 SNu (Tu- 
ratto, Cappellaro, & Benetti 1994), where h = H o /10Q 
and 1 SNu = 1 SNIa per century per lO 10 ^©- 

In principle, constraints on the SNIa rate can be 
given also by estimates of the iron abundance in galac- 
tic flows (see, e.g., Rcnzini et al. 1993, and references 



therein). These were first attempted using data from 
the Ginga satellite for NGC4472, NGC4636, NGC1399 
(Ohashi et al. 1990; Awaki et al. 1991; Ikebe et al. 1992), 
and then from the BBXRT sa tellite for NGC1399, 
NGC4472 ( |Serlemitsos et al. 1993| ), and more recently 
from ASCA with a superior spectral energy resolution 
( Locwcnstcin et al. 1994j ; Awaki et al. 1994; Arimoto et 
al. 1997; Matsumoto et al. 1997). Under the assumption of 
solar abundance ratio, the analysis of all these data sug- 
gests a very low iron abundance, consistent with no SNIa's 
enrichment and even lower than that of the stellar compo- 
nent. However, some authors have found that more com- 
plex multi-temperature models with higher abundance 
give a better fit of the data ( Kim et al. 1996 , Buote & 
Fabian 1997). Moreover, the above results are based on 
iron line diag nostic tools whose r eliability has been ques- 
tioned (e.g., Arimoto et al. 1997), especially because of 
uncertainties affecting the Fe L-shell atomic physics for 
temperatures less than 2 keV, that are typical of hot gas 
flows in ellipticals (Liedahl, Osterheld, & Goldstein 1995). 
Line emission even in simple (e.g., isothermal) astrophys- 
ical plasmas needs to be understood, and reliable fits to 
the data made, before secure consequences concerning the 
SNIa rate can be drawn from X-ray determined abun- 
dances. 



In summary, recent optical studies agree on a present 
epoch SNIa rate much lower than previously used, and in- 
dicate that also the dark matter content could be lower. 
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Moreover, cuspier density profiles, especially for the stellar 
distribution, should be used. All this raises the question of 
whether the CDPR scenario is altered by these changes in 
the main ingredients of the problem, and more in general 
what are the effects on the properties of hot gas flows. 
In this paper we address these points, with a new set of 
hydrodynamical simulations. In Section 2 we present the 
galaxy models, the source terms, and the integration tech- 
niques. In Section 3 we discuss the main properties of the 
evolution of gas flows in our new models, and we compare 
them with the observations and the CDPR results. In Sec- 
tion 4 we discuss the results using energetic arguments, 
and in Section 5 the main conclusions are summarized. 



2. The models 

2.1. Galaxy models 

Following the discussion in Sections 1.1 and 1.2, our model 
galaxies are a superposition of a Jaffe density distribution 
for the luminous matter p* [eq. (1) with M = M„, r c = r*, 
7 = 2], and of a Hcrnquist distribution for the dark mat- 
ter [eq. (1) with M = M h , r c = r h , 7 = 1]. We call 
these two-component mass models JH models. Both den- 
sity distributions have finite mass, and so no truncation 
radius is applied. The relation between the scale length r* 
and the effective radius for the Jaffe model is R c ~ 0.76r*. 
We constrain the luminous bod y to lie on the fundamen- 
tal plane of elliptical galaxies ( Djorgovski fc Davis 1987 ; 
Bender, Burstein, & Faber 1992). The input observables 
are Lb, and the central value of the projected velocity dis- 
persion o~ c . The stellar mass is related to o~ c through the 
projected virial theorem, as described in Ciotti, Lanzoni, 
& Renzini (1996); the luminous matter distribution is then 
completely defined, because it is assumed that the stellar 
mass-to-light ratio is constant with radius, and the orbital 
distribution is isotropic. The dark matter distribution is 
determined by choosing the ratios of the dark to luminous 
mass (1Z — Mh/M*), and of the dark to luminous core 
radius {[3 — r^fr*). The basic dynamical properties of the 
JH models (radial trend of the velocity dispersion, kinetic 
and potential energy, etc.) are given in the Appendix. 

2.2. Source terms 

The time evolving input ingredients of the numerical sim- 
ulations are the rates of stellar mass loss, of SNIa heat- 
ing, and of thermalization of the stellar velocity dis- 
persion; these are calculated as in CDPR, where a de- 
tailed description is given. Here we summarize only the 
main properties of the input quantities. The stellar mass 
loss rate is accurately described by M*(i) ~ 1.5 x 
10~ n Lb t\s M&yr~ l , where Lb is in L Q , and tis is 
time in units of 15 Gyr; in the numerical code the exact 
mass return prescribed by the stellar evolution theory is 
used. The SNIa heating rate is parameterized as Lsn(<) = 
%%(i)£B = 7.1 x IO^nLe^ 1 ' 5 er S s_1 > where 



Fig. 2. The cumulative mass profiles normalized to the total 
stellar mass M„ of JH models with 1Z — 2 and 1Z = 9; the cases 
(3 = 2 and (3 = 4 are shown in each panel. Within r — 2R C , 
1Z = 0.6 and 0.3 respectively if /3 = 2 and j3 = 4 (upper panel); 
11 = 2.8 and 1.1 if f3 = 2 and (3 = 4 (lower panel). 



Lsn = 10 51 erg is the kinetic energy injected in the ISM 
by one SNIa, Rm{t) = i9 SN 0.88 h 2 i^ 1 ' 5 SNu, and h = 0.5. 
When $sn = 1 and t\s = 1 i?sN is the rate estimated by 
Tammann (1982). LgN is assumed to be decreasing with 
time just faster than the stellar mass loss rate. Unfortu- 
nately, given our enduring ignorance of the nature of the 
SNIa's progenitors (e.g., Branch et al. 1995) the evolution 
of the SNIa rate is not known; some arguments though, as 
an observed pronounced correlation of the rate with the 
star formation rate along the Hubble sequence, favour a 
declining rate (see also Renzini et al. 1993; Ruiz-Lapuente, 
Burkert, & Canal 1995; Renzini 1996). With the assump- 
tions above the specific heating for the gas is decreasing 
with time, and so we have the wind/outflow/inflow sec- 
ular evolution in the CDPR models (§1). This direction 
of the evolution of the flow is maintained as long as Rsn 
decreases faster than AL, and the results do not quali- 
tatively depend on the particular slope adopted for Rsn', 
when this slope is much higher than that of M*, though, 
too much iron is produced at early times (see Renzini et 
al. 1993 for a more extensive discussion). 

The heating given by the thermalization of the stellar 
velocity dispersion at each radius is determined using the 
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velocity dispersion profile obtained by solving the Jeans 
equation in the global isotropic case (see the Appendix). 

2.3. Numerical simulations 

The time-dependent equations of hydrodynamics with 
source terms, and the numerical code used to solve them, 
are fully described in CDPR, together with the initial con- 
ditions. We adopt here a higher spatial resolution in the 
central regions of the galaxies: the central grid spacing is 
30 pc instead of 100 pc. This allows a better sampling of 
the inner regions now described by a power law, and is ad- 
equate to highlight the differences in the flow properties 
between JH models and King models plus a quasi isother- 
mal dark halo. The model galaxies are initially devoid of 
gas, a situation produced by the galactic winds established 
by type II supernovae, early in the evolution of elliptical 
galaxies. The gas flow evolution is followed for 15 Gyrs. 



Fig. 3. The evolution of the number density n, the temper- 
ature T and the velocity v of the gas, in a typical PW case 
(the model with Lb = lO n L , TZ = 2, (3 = 4 and ??sn = 0.25 
in Table 1). The solid line refers to t = 1.5 Gyrs, the dotted 
line to t — 4.5 Gyrs, and the dashed line to t = 15 Gyrs. In 
this model the initially supersonic inflow slows until it becomes 
fully subsonic at t ~ 6 Gyrs; the outer sonic radius starts at 
~ 2R e and moves outward, until it establishes at ~ 8R e from 
t ~ 10 Gyrs on. 



3. Evolution of the gas flows 

We present here the results of the hydrodynamical simu- 
lations. The basic input parameters and output quantities 
are shown in Table 1 . Two representative blue luminosities 
have been used to investigate the typical gas flow behavior, 
Lb = 5 x 1O 1O £0, and Lb — 10 11 L Q ; the corresponding 
M* are 1.9 x lO n M and 4.2 x lO n M . Various values of 
(TZ, (3, #sn) are chosen. #sn = 0.25 is close to the most se- 
cure current estimate coming from optical surveys (which 
is i?sn ~ 0.27, see §1.3); to consider the indications com- 
ing from the low iron abundances given by X-ray data 
(§1.2), we also explore the case -#sn =0.1; for comparison 
with the long-lived inflow case, we also use i?sn = 0. The 
distribution of the dark mass is broader than that of the 
luminous mass (/3 = 2 or = 4); the dark mass prevails 
outside one or two R c if TZ = 9, or outside a few R e if 
TZ = 2 (Fig. 2). 

A general feature of all the evolutionary sequences, ex- 
cept for those with $sn = 0, is to develop a central inflow 
at early times, while the outer regions are still outgassing, 
i.e., a decoupled gas flow (see Fig. 3). The inflow region is 
bordered by a stagnation radius r st that may range from 
a small fraction of R e to many R e , at the end of the evolu- 
tion (see Table 1). If r s t maintains its position within a few 
R c , we call the flow partial wind (hereafter PW; see the 
representative case in Fig. 3 and 4); if it increases con- 
siderably to more than 10i? e , we call the resulting flow 
global inflow, even though this is strictly the case only 
when #sn = 0. In the PW case, radiative losses suppress 
an outflow in the inner parts of the galaxy, causing a small 
inflow region, but a wind can be sustained in the external 
parts, where the gas density is much lower and the gas is 
also less tightly bound. The X-ray luminosity of PWs and 
inflows steadly decreases with time, following the decrease 
of M* and of the heating. Note that the stagnation points 
resulting from the simulations are well outside the first 
grid point, so the flow is properly resolved. 

Another general trend in the results is that, at fixed 
TZ, as $sn and/or (3 decrease, the larger are the inflow 
region, Lx, the mass accreted to the center M- m , the mass 
flowing to the center per unit time M- m , and the smaller is 
M out , the mass escaping the galaxy per unit time, because 
the heating is lower and/or the dark mass concentration 
is higher. 



3.1. L B = 5x 10 10 L© models 

For TZ = 2, and $sn > 0, these galaxies develop PWs with 
r st varying from 0.02 R c to ~ 3i? . When $sn > 0.25 a 
global wind persists until the present time, and so Lx is 
very low. Even when z?sn = 0.25 the hot gas, in PW, shows 
such a low Lx that the observed X-ray emission must be 
largely produced by the stellar sources, whose emission has 
been estimated to be of the order of ~ 10 40 erg s _1 (see 
Fig. 1). When t? S N = 0.1, the final L x < 5 x 10 40 erg s"\ 
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some of them could be explained with accretion effects 



(see, e.g., 


Rcnzini ct al. 1993 


Kim & Fabbiano 1995). 



Bcrtin & Toniazzo 1995 



Note that for PWs (TZ = 2) L x is higher when i? SN 
is lower, because the inflow region is larger. For global 
inflows (TZ = 9) Lx is higher when $sn is also higher, 
because the gas is hotter. Altogether, the variation in Lx 
in Table 1 is a factor of ~ 100; it is of a factor of ~ 30 
when TZ = 2, and just ~ 2 when TZ = 9. 



Fig. 4. Upper panel: the time evolution of the mass loss M„, 
of the mass inflowing into the central sink Mi n , and of the 
mass lost by the galaxy M out . Lower panel: the time evolution 
of the stagnation radius r s t. All refers to the representative 
model chosen for Fig. 3. 



and the hot gas starts to dominate the X-ray emission of 
the galaxy. Significantly higher Lx values are obtained by 
increasing the gas temperature, not just by retaining more 
gas; in fact the model with TZ = 2 and $sn = 0, a global 
inflow that lasts for the whole galaxy lifetime, has just 
Lx ~ 7 x 10 40 erg s _1 . To increase the gas temperature 
and Lx, the dark matter content must be increased first, 
because global winds are obtained by just increasing the 
SNIa heating. A deeper potential well by itself produces 
hotter gas, because of a higher gravitational compression 
during the inflow; moreover, it can be coupled to a larger 
inflow region, if i9§n is kept constant, or to a higher heating 
without full degassing if i9gx is increased. 

For TZ = 9 and $sn^0.25 the gas is always in a 
global inflow at the present time. The highest Lx ~ 
2 x 10 41 erg s _1 is obtained when ^sn = 0.25 and = 2, 
i.e., when the heating by SNIa's and gravitational com- 
pression is the highest possible for the chosen ranges of 
these parameters. This Lx is still lower than the high- 
est Lx observed at L B = 5 x 10 10 L Q (Fig. 1), a prob- 
lem that was known also from previous works (see the 
references cited in the Introduction). However, Lx as 
high as ~ 10 42 erg s _1 are comparable to those typi- 
cal of poor clusters rather than single galaxies; at least 



3.2. L B = 10 n L Q models 

Models with TZ = 2 again develop PWs, with quite larger 
central inflow regions than above for the same range of 
$SN and (3 values, due to a deeper potential well. The 
latter makes also the gas hotter. Larger inflow regions 
and higher gas temperatures, coupled to a larger mass 
return rate, make the X-ray emission in these galaxies al- 
ways higher than that of the models described in §3.1, for 
the same (TZ, (3, i?sn)- Lx goes from 3 x 10 40 erg s _1 to 
2 x 10 41 erg s _1 , when the flow is a PW. So, it is always 
comparable to or higher than the estimate of the stellar 
X-ray emission at Lb = 10 11 L Q (Fig. 1); this is in agree- 
ment with the fact that in Fig. 1 all galaxies lie above this 
estimate, for Lb^W^Lq. 

As expected, we always have global inflows ii TZ = 9, 
and the highest Lx is reached when TZ = 9 and $sn = 
0.25, i.e., again the most efficient way of increasing Lx 
when the flow is not decoupled is to increase TZ, and then 
$sn- The spread in Lx is of a factor of ~ 10 if TZ — 2, and 
very small for global inflows. The total spread is roughly 
of a factor of 30, lower than for Lb = 5 x 1O 1o L . 

3.3. Other models 

In order to investigate how common is the PW phase, 
we have run two more series of evolutionary sequences 
corresponding to a low Lb = 2.2 x 10 10 L Q , and to 
Lb = 7.1 x 10 10 L Q ; moreover, also the case TZ = 4 has 
been explored. The results are summarized in Table 2. 
PWs populate again most of the parameter space. The X- 
ray luminosity of Lb = 2.2 x 10 10 L Q models remains at 
low values even at high TZ, in agreement with the obser- 
vations (Fig. 1); for many of these galaxies Lx is likely to 
be dominated by the stellar emission. 

3.4- Emission temperatures and X-ray surface brightness 
profiles 

Emission temperatures of the hot gas have been calculated 
recently by Matsumoto et al. (1997) and Buote & Fabian 
(1997), using ASCA data and two-temperature model fit- 
ting, for about 20 early-type galaxies. These temperatures 
cover the range ~ 0.3 — 1.0 keV, and the bulk of values 
lies within 0.5 — 0.8 keV. The X-ray luminosity weighted 
emission temperatures of the models are given in Tables 1 
and 2. They lie in the range ~ 0.3 — 0.8 keV, and so are in 
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Table 1. Results of the JH model evolution. Lx is calculated for the (0.2-4) keV band; T is the emission weighted temperature 
of the flow. Mi n is the gas mass flown to the center in 15 Gyrs; Mm and M ou t are the gas masses that flow to the central sink 
and are lost from the galaxy, per unit time, at 15 Gyrs. At this time M* = 0.72 M Q /yr, for the model with Lb = 5 x 10 10 L©, 
and M» = 1.44M Q /yr for the model with L B = 10 11 i Q . \ is defined in §4. 



good agreement with those observed. They increase with 
central velocity dispersion, dark matter mass, and SNIa 
rate. Note how a substantial SNIa's heating helps obtain- 
ing temperatures close to those observed; other possible 
sources of heating not studied here are the presence of 
an external pressure due to an intracluster or intragroup 
medium (e.g., Bertin & Toniazzo 1995), and that of a cen- 
tral black hole (Ciotti & Ostriker 1997). 

The comparison with the observed X-ray surface 
brightness profiles Sx(-R) is more delicate, because there 
is not a typical profile, nor a range of typical profiles. In 
the next, first we examine what kind of Sx(-R) are shown 
by PWs (which are the resulting flow regimes of the bulk 
of the galaxies, in our scenario), and then we discuss how 
Ex(-R) of our models compare with the available observa- 
tions. 

In Fig. 5 we compare the shapes of the X-ray profiles 
of partial wind solutions with different stagnation radii 
with that of a global inflow with $sn = (which closely 
resembles the cooling flow solution). When the stagnation 
radius is larger than the optical effective radius, the shape 
of Ex (R) is indistinguishable from that of the inflow, over 
most of the galaxy. PWs with very small stagnation radii 
are instead less steep outside ~ R c ; this is because in the 
external regions the gas density is actually lower for the 



PW than for the inflow, but the prevailing effect is given 
by the higher gas temperature of the PW, produced by 
the higher $sn- So, judging from the shape of the Sx(-R) 
profile, inflows and PWs should be indistinguishable by 
X-ray observations, if r st is a few R c . Of course, when r st 
is very small, Lx is very low too, and so these PWs are 
much less detectable by X-ray observations than global 
inflows. 

Now let's turn to the comparison with the observa- 
tions. When Lx of the gas flow is very low, and the X-ray 
emission is dominated by the contribution of the stellar 
sources, our models exhibit the X-ray profile of the un- 
derlying stellar population. No detailed Ex(-R) for low 
Lx/Lb galaxies is available yet from the observations. 
When Lx of the flow is high, the observations - espe- 
cially those produced by ROSAT - show a large variety 
of shapes for T,x(R), often modified also by the interac- 
tion with the environment (e.g., Trinchieri, Fabbiano & 
Kim 1997). Using Einstein data, Trinchieri, Fabbiano, & 
Canizares (1986) found that the observed S x (i?) tend to 
follow the optical ones, for few best studied X-ray bright 
galaxies. In our scenario such galaxies host PWs with large 
stagnation radii, or inflows; these flow phases have the 
same problems as cooling flows to reproduce X-ray pro- 
files following the optical ones: they are too peaked (see 



S. Pellegrini and L. Ciotti: Decoupled hot gas flows in elliptical galaxies 



8 



Table 2. Results of additional JH model evolutions. All quantities are defined as in Table 1. After 15 Gyrs M 4 = 0.32 Mq/yt, 
for the model with L B = 2.2 x 10 10 L Q , and M, = 1.02 M Q /yr for the model with L B = 7.1 x 10 10 L©. At fixed L B , the models 
with the larger R c are less concentrated than required for them to lay exactly on the fundamental plane. 



Fig. 5). They can be brought in agreement with the ob- 
served Sx(-R) by the introduction of distributed mass de- 
position due to thermal instabilities in the hot g as, clS 
shown and discussed in detail by Sarazin fc Ashe 1989] 
and Bertin & Toniazzo 1995; another possibility is sta- 



tiona ry convcctive accretion onto a central massive black 
hole ( Tabor fc Binney 1993 ), or unstable accretion (Ciotti 
& Ostriker 1997). 

A detailed comparison with Y.x(R) observed for a spe- 
cific galaxy requires a large number of simulations to find 
whether there is a combination of input parameters (with 
the possible addition of the effect of the environment, 
and/or of a central black hole) that gives a model re- 
producing the observations. This goes beyond the scope 
of this paper; it has been done successfully for NGC4365, 
in the framework of the CDPR scenario, by Pellegrini & 
Fabbiano (1994). 

3. 5. Effects of the new ingredients 

We summarize here the main properties of the flows, and 
the main differences with the results obtained by CDPR, 
produced by the new mass distributions and by the re- 
duction of i9sn and 1Z. The global wind phase disappears 
for all galaxies except the smallest ones (see the case 
L B = 2.2 x 10 10 i© in Table 2): a central inflow, even 
though very small, is always present from the beginning are needed (Awaki et al. 1994); ASCA data require cen 



of the evolution. As a consequence the outflow phase - the 
transition from a global wind to a global inflow in CDPR 
- also disappears. Large Lx variations are produced in 
the JH models by the different size of the central inflow 
region. Most of the observed spread in Lx at fixed Lb can 
be reproduced again, as it was in CDPR (Fig. 1): besides 
the observed variation from galaxy to galaxy in o~ c and 
in the concentration R e , a spread in the dark matter con- 
tent, and/or a variation in the SNIa rate can produce large 
variations in r st and in the gas temperature, and then in 
L x (of even a factor of - 100, Table 1 and Fig. 1). If ^ 
remains at low values (TZ<4), though, it is difficult to re- 
produce Lx higher than a few times 10 41 erg s _1 (Table 
2), and so to cover all the observed Lx variation. More- 
over, for a given range of variation of 1Z and i?sn, more 
luminous galaxies show less scatter in Lx- 

Another property of the flows in JH models is to accu- 
mulate cold gas at their centers, during the evolution (M; n 
in Tables 1 and 2). This mass when #sn > ranges from 
10 9 M Q to 2 x 10 11 M Q . Evidence of cold gas at the center 
of bright X-ray galaxies has come recently from the high 
column densities (Nh) required to obtain good fits of their 
X-ray spectra. The B BXRT data require ~ 1 10 M Q of 
cold gas in NGC1399 ( |Scrlcmitsos ct al. 1993]) ; to fit the 
ASCA data of NGC4472, 3 x 10 10 M m of cold material 
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potential well (L g ), and the heating supplied by the SNIa 
explosions plus the thermalization of the stellar velocity 
dispersion (L a ). x was introduced by CDPR to predict the 
expected flow phases in King galaxies, and was shown to 
work remarkably well by the numerical simulations: when 
X > 1 the flow turned out to be an inflow, when x < 1 a 
wind. In Tables 1 and 2 most of the resulting flow phases 
are not global, but PWs, and so we are moved to investi- 
gate in more detail the distribution of the energetics inside 
the galaxy models. This distribution in fact can be much 
different for two galaxies with the same global x- The local 
Xe is defined as 

dLJr) 

Xi{r) = „ - f ,l' J T — f-r, 2 
where 

dL g (r) = 47ra(i)p*(r)|0 T (r)|r 2 dr, (3) 



Fig. 5. The X-ray surface brightness profile of three models 
with L B = 5 x 10 10 L Q , K = 2 and (3 = 2: the solid line 
corresponds to a global inflow ($sn = 0), the dotted line to 
a PW with r«t = 3.5i? c (the model with i?sn = 0.1 in Table 
1), and the dashed line to a PW with r st = 0.04R C (the model 
with $sn = 0.25 in Table 1). The heavy dot-dashed line is the 
optical surface brightness profile in arbitrary units. 



tral column densities larger than the Galactic value also 
for NGC46 36, NGC4406, NGC7 20, NGC1399, NGC1404, 
NGC4374 (|Arimoto et al. 1997[). It has been suggested to 



explain this finding as the accumulation of cold material 
by a steady state cooling flow. This accumulation could 
have been produced by a PW as well; moreover, even dif- 
ferent amounts of cold material can be naturally explained 
in this scenario (see the large variations of M- m in Table 1). 
So, the finding of cold matter at the center of early- type 
galaxies cannot be considered as evidence for a long lived 
cooling flow. Note however that the accuracy of these Njj 
absorption measurements could be affected by the uncer- 
tainties plaguing plasma spectra, and producing the iron-L 
problem (§1.3). 

4. An energetic explanation for the occurrence of 
the PW in different galaxy models 

The origin of PWs, and of the flows of a different na- 
ture obtained by CDPR, can be explained qualitatively in 
terms of the local Xi function. The global x is defined as 
the ratio between the energy required to steadily extract 
the gas lost by the stars per unit time from the galactic 



dL a {r) = 6ira(t) p^{r)a1{r)r 2 dr, 



and 

d L SN (r) = ^^4^(r)r 2 dr, 



(4) 



(5) 



where 4>t = <t>* + 4>h is the total potential, o~„ is the one- 
dimensional isotropic velocity dispersion obtained by solv- 
ing the Jeans equations fo r the adopted mass model (e.g., 
Binncy fc Trcmainc 1987 ), and a(t) = M*/M* is the spe- 
cific mass return rate. The global x is obtained by replac- 
ing the differentials in equation (2) with their integrals 
over the whole galaxy. 

The differences in the flow behavior of JH models and 
models with a central constant density region can be un- 
derstood with the aid of Fig. 6, where xi is shown for 
four representative one-component models. These have the 
same global x, and so from an energetical point of view 
are globally equivalent. Three of them are 7-models, with 
7 = 2,1,0, i.e., they are the Jaffe model, the Hernquist 
model, and a model with a core. The fourth density dis- 
tribution is described in the Appendix, and is flatter at the 
center than the 7 = model, but for large r it decreases 
oc r~ 4 , as all the 7-models. This distribution is studied 
for its similarity with the King (1972) model, that needs 
the annoying introduction of a truncation radius, and has 
more complicated dynamical properties. 

It is apparent from Fig. 6 that the steeper the density 
profile, the stronger is the variation of xe across the galaxy, 
and the higher are the values that it reaches at the center. 
This explains why a strong decoupling can be present in 
the flow of highly concentrated systems: a xi significantly 
higher than unity in the central regions produces a central 
inflow, while in the external parts a degassing is energet- 
ically possible because \t < !• This can happen when the 
global x is either > 1 or < 1 (see Fig. 6 and Table 1), so 
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order of the core radius of the dark halo), while with low 
1Z this region is of the order of the stellar core radius. 

In summary, the combined effect of SNIa's and dark 
matter is more varied for models with a core than for JH 
models: for a large range of 1Z values, the latter keep in 
PWs with a varying r st , while the former can be in wind, 
PWs, outflows or inflows. 



Fig. 6. The radial trend of %t as a function of r/R e for the 
models considered in §4. The ratio R c /r* is given by 0.76, 
1.82, 2.87, 1.73 respectively for the Jaffe (solid line), Hernquist 
(dotted), 7 = (short-dashed) and the model described in the 
Appendix (long-dashed). 



this parameter is not a good indicator of the flow phase 
for highly concentrated systems. If x decreases, the radius 
where xi = 1 moves inward for all the models, and so a 
larger part of the galaxy is degassing. Note a fundamental 
difference in the trend of Xi( r ) as a function of x for mod- 
els with a core and cuspier models: while in cuspy models 
there is always a region with xt > 1) even though small, 
this is not the case for core models, in which such a region 
suddenly appears by increasing with a size larger than 
the core. 

A second effect on xi is produced by the dark halo. If 
this halo is more diffuse than the stellar mass, increasing 
7Z makes more bound the external regions, when x is kept 
constant, and so the radius at which Xi = 1 moves outward 
(see Fig. 7 for the case of JH models). 

The trend of Xi{ r ) f° r models with a core explains 
the results of the numerical simulations for King mod- 
els plus diffuse quasi-isothermal dark halos obtained by 
CDPR and Pellegrini & Fabbiano (1994). The flow phases 
found by CDPR with 1Z = 9 were all global; with low 
7Z, instead, the flow can be decoupled again, as found by 
Pellegrini & Fabbiano (1994) in their detailed modeling of 
the X-ray properties of two ellipticals. This because with 
high 1Z the region with xe ~ constant is very large (of the 



Fig. 7. The radial trend of \t as a function of r / R c for the JH 
model. The different curves are labelled with the value of 1Z, 
while the global \ an d P are constant at the values indicated. 



5. Conclusions 

Using evolutionary hydrodynamical simulations, we have 
investigated the properties of hot gas flows in elliptical 
galaxies described by a Jaffe stellar density profile, and 
containing various amounts of dark matter, distributed 
as a Hernquist law; the rate of SNIa explosions is one- 
fourth of the Tammann's rate, or lower. The results can 
be summarized as follows. 

The simulations show that most of the gas flows are 
strongly decoupled, i.e., they develop an inflow in the cen- 
tral region of the galaxy, while the external parts are still 
degassingt!l. The stagnation radius of these flows (PWs) 

Note that decoupled flows are present also in flat and ro- 
tating models with low i?sn, albeit for reasons different from 
those causing the PWs studied here (D'Ercole & Ciotti 1997). 
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may range from a small fraction of the optical effective 
radius to many R e . This situation is present even when 
the global energy balance of the flow indicates that the 
energy inputs from SNIa's and from the thermalization of 
the stellar random motions is high enough to unbind all 
the gas (x < 1)- Alternatively the gas can be outflowing 
from the outer part of the galaxy even when the global 
energy balance indicates that the energy available is less 
than needed to unbind it all (x > 1)- So, \ 1S n °t a good 
indicator of the flow phase; larger x values, though, cor- 
respond to larger stagnation radii, for a fixed Lb- Global 
inflows are produced when the dark matter content is so 
high to bind the gas over the whole galaxy. 

The X-ray luminosity of PWs is higher when $sn is 
lower, because the inflow region is larger; for global inflows 
L x is higher when -d SN is also higher, because the gas is 
hotter. The key factor causing large Lx variations in JH 
models is the size of the central inflow region. The lowest 
Lx/Lb values observed can be easily reproduced by PWs; 
high dark matter contents are required to approach the 
highest Lx observed. 

The shape of the X-ray surface brightness profile of 
PWs is close to that of global inflows, over most of the 
optical image, if the stagnation radius is a few R c ; it is 
externally less peaked than that, if r st is a fraction of R c , 
i.e., when the gas is hotter than in the inflow case. 

The strong decoupling of the flow is explained by the 
radial trend of the local energy balance xti that is very 
peaked for steep mass density profiles as in the JH mod- 
els. Previously used density profiles which were flat at the 
center - such as the King models plus quasi-isothermal 
dark halos - showed a larger tendency to have global flow 
phases, because also their xt( r ) is flatter. So, the change 
of the mass profile has an important effect: both numerical 
simulations and analytical calculations show how peaked 
density profiles preferentially develop decoupled gas flows 
(in the sense that the region of the parameter space pop- 
ulated by PWs is large compared to that of inflows or 
winds). 

The results of our simulations compare with the obser- 
vations in two main aspects. The first is that most of the 
observed spread in Lx shown by the data is easily repro- 
duced: the highest Lx are again, as in CDPR, associated 
with global inflows, while the bulk of the galaxies are now 
in PW, rather than in outflow; global winds are present 
only in the smallest galaxies. 

The second is that the presence of cold gas at the 
center of Es, and a peaked X-ray surface brightness pro- 
file, cannot be unequivocally associated to a cooling flow: 
a PW could be present as well, with a significant part 
of the galaxy degassing. Particularly, a galaxy can have 
$SN = 0.25 and still a non negligible amount of cold mat- 
ter at the center. In addition, since r st varies over a wide 
range of radii, PWs can accumulate largely varying quan- 
tities of cold gas. 
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6. Appendix 

6.1. The J H models 

We derive here the main properties of the JH models 
that have been used in the numerical simulations and 
in the computations of x and Xi- The reported quanti- 
ties are dimcnsionlcss, and the normalization constants 
are expressed in terms of the total stellar mass M», the 
core radius r* and the gravitational constant G; again 
1Z = Mh/M* and (3 = r^fr*. The cumulative mass inside 
s — r j r ^ anc i total potential are given by: 



M(a) 



+ 



TZs 2 



l + s (0 + s) 



2 ' 



<M S ) = - In 



K 



(6) 



(7) 



In the hydrodynamical simulations we used the globally 
isotropic stellar velocity dispersion profile, that we cal- 
culate here. For any two-component model the velocity 
dispersion profile is given by p*cr 2 = J**(s) + TZI^(s), 
where 0% is the 1-dimensional radial velocity dispersion. 
The first term is the isotropic velocity dispersion for the 
Jaffc model: 



-^** (5) — 



(6s 2 + 6s- l)(2s+ 1) 
8tts 2 (1 + s) 2 



A ln I± 

2tt s 



(8) 



and the second one describes the contribution due to the 
Hernquist dark halo: 



+ 



lm»Qg + l) _ ln(l + s)Q8-2) _ lnQ3 + s)(2/? - 1) 

2tt/? 3 2tt(/?- 1) 3 2tt/? 3 (/3- l) 3 

2Q3 2 -0 + l)s 2 + {[3 + 1)(2/? 2 - 3/3 + 2)s + /3(j3 - l) 2 
Airf3 2 {(3 + s){l + s) 



(9) 



Physically acceptable limits can be easily found for (3 = 
(thus mimicking a black hole at the center of a Jaffc 
model) and for (3 = 1. 

The cumulative potential and kinetic energies inside r 
can be found analytically for this two-component model, 
but for simplicity we report here only their total values. 
The total potential energy is given by J7* = ?7** + 1ZU*h, 
where 



U** - 
and 
U*h - 



1 

2' 



(/3-1) 2 ■ 



(10) 



(11) 
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By integrating over the whole galaxy equation (3), we 
have L g = a(i)(GM 2 /r*)(2|i7**| + 1l\U* h \). The kinetic 
energy is again the sum of two different contributions, 
= K** + lZK*h, where 



and 



1 

4' 



(/3 + l)ln(/3) + 2(1-/3) 



2(/?-l) 3 



(12) 



(13) 



The integration over the whole galaxy of equation (4) gives 
La = a(t)(GM^/r,){K„+TZK, h ). 

6. 2. A centrally flat model 

Here we derive the dynamical quantities corresponding to 
the density distribution 



P*( r ) = 



1 



1 



7T 2 (1 + S 2 ) 2 ' 

The mass contained inside r is given by 



M*{s) 



arctan(s) 



1 



(14) 



(15) 



and the potential is 
- 2arct^l_ 

TT S 

The globally isotropic velocity dispersion is given by: 

3arctan(s) 2 (2 + 3s 2 ) arctan(s) 4 + 3s 2 



2tt 3 



s(l + s 2 )ir 3 



2(1 + S 2 ) 2 7T 2 



As for the JH model, the potential and the kinetic energy 
inside the radius r can be expressed analytically, but here 
we give only their total values, being the only ingredient 
required in the definition of x- 



2tt' 



(18) 



and 



l 

47r' 



(19) 



Finally, the projected distribution associated to is: 
1 



2tt(1 + j R 2 ) 3 / 2 ' 
from which one obtains R e = V%r* 



(20) 
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